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Master equation for a quantum particle in a gas 
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The equation for the quantum motion of a Brownian particle in a gaseous environment is derived 
by means of S-matrix theory. This quantum version of the linear Boltzmann equation accounts non- 
perturbatively for the quantum effects of the scattering dynamics and describes decoherence and 
dissipation in a unified framework. As a completely positive master equation it incorporates both 
the known equation for an infinitely massive Brownian particle and the classical linear Boltzmann 
equation as limiting cases. 
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How is the motion of a quantum particle affected 
by collisions with ambient gas molecules? The well- 
established field of Quantum Brownian Motion 0,0] pro- 
vides an answer provided the interaction can be linearized 
and the particle state is close to classical. One is faced 
with a rather different situation if the Brownian particle 
is in a highly non-classical motional state, say, due to 
the passage through an interferometer or the entangle- 
ment with another controlled degree of freedom. Such 
correlations are becoming experimentally accessible in 
the emerging field of molecular quantum optics, where 
the quantum nature of molecular motion is tested and 
exploited 3]. 

In order to assess the (partial) loss of coherence in 
the case of strong, nonclassical correlations in the mo- 
tional state it is necessary to provide a detailed, non- 
perturbative account of the microscopic scattering pro- 
cess. At the same time, for the effects of decoherence 
to be relevant at all, the gas density must be sufficiently 
low, so that the environmental gas can be safely taken 
as not self- interacting and Markovian Q . In analogy to 
the classical case the associated description of a sin- 
gle distinguished tracer particle within an ideal gas may 
be called a linear Boltzmann equation. It should not be 
confused with the linearized Boltzmann equations for the 
reduced single particle gas state, obtained perturbatively 
from its multi-particle description. 

The investigation of the loss of coherence due to gas 
collisions was initiated by Joos and Zeh who con- 
sidered the limiting case of an infinitely massive tracer 
particle. This theory, which was later refined 0, HI and 
tested experimentally describes the pure spatial "lo- 
calization" of an extended coherent matter wave into a 
mixture with reduced spatial coherence, but it cannot 
account for dissipation. The situation is much more in- 
volved if the tracer mass M is finite and comparable to 
the gas mass m so that the ratio m/M must not be ne- 
glected. In this case "localization" occurs both in posi- 
tion and in momentum |lCll |. and the appropriate kinetic 
equation must describe the full interplay of decohering 
and thermalizing dynamics. 

So far, the most important advancement in this direc- 



tion is the proposal by Diosi of an equation based 
on a combination of scattering theory and heuristic ar- 
guments. A more recent development is the theory by 
Vacchini in terms of the dynamic structure factor of 
the medium, an approach limited to the (weak coupling) 
Born approximation, like those in 0, Q|. 

This letter presents the full quantum version of the 
linear Boltzmann equation, describing the whole range 
of coUisional effects from decoherence to dissipation. It 
provides a transparent and stringent derivation, and a 
discussion of its implications and limits. The only essen- 
tial premise is the Markov assumption, which implies — in 
the spirit of Boltzmann's classic derivation — that both 
the rate and the effect of individual two-particle scat- 
tering events are separately physically meaningful, while 
subsequent collisions with the "same" gas molecule are 
negligibly unlikely. Accordingly, the gas may be taken 
to be ideal (not self-interacting), stationary (diagonal in 
momentum), and uniform in position space, thus cover- 
ing e.g. thermal states of bosons and fermions, but no 
liquids. 

In operator form the master equation reads dtp — 
{ih)~^ [H,p] + Cp with H = PV(2M) the free Hamil- 
tonian and 



Cp = 



-QK 



-2P^Q,K^Q,K - ^^q,K^Q,KP 



(1) 



Here the integration is over all momentum transfers Q, 
and for fixed Q also over the perpendicular plane = 
{K e ]R3 : X • Q = 0}. The Lindblad (jump) operators 
have the form 



-Q,K 



e^R-Q/^F (fC, P; g) 



(2) 



where R and P are the position and the momentum op- 
erator of the Brownian tracer particle. The function F, 
which is operator- valued in contains all the details 
of the coUisional interaction with the gas. It involves the 
elastic scattering amplitude / (Pout j EiaJ > the momentum 
distribution function /i (p) of the gas|2flj, and its number 



2 



density rigas- It is convenient to denote relative momenta 
by 



rel p,P) p —P, 

m M 



with m* = mM/ {M + m) the reduced mass. More- 
over, for given Q ^ let us denote the parallel and the 
perpendicular contribution of a vector (operator) P by 
P\\Q = {P-Q) Q/Q^ and by P^q = P P||q, respec- 
tively. With these definitions 



F{K,P;Q) 



(3) 



X/ ( rel {K^Q, P^q) - |, rel (fsT^g, P^q) + | 



1/2 



This implies that both the scattering amplitude and the 
distribution function attain an operator character in ij^J, 
and that the particle momentum contributes with the 
part perpendicular to the momentum exchange to the 
former and with the parallel one to the latter. For phys- 
ically reasonable interactions the function F decreases 
sufficiently fast as Q ^ so that |^ is well-defined. 

Note that the form of £p fits the general structure 
of a translation-invariant and completely positive master 
equation, as characterized by Holevo [l^ (se e llal for a 
discussion), although the summation in Ref. '1^ is here 
replaced by the integrations in We will see below 
that the master equation assumes a more intuitive form 
in the momentum representation. 

My first aim is to provide a derivation of Q)-©. To 
that end, let us first define the positive operator f which 
yields the total collision rate. As in classical mechanics, 
the rate is determined by the gas density, the modulus 
of the relative velocity i; = |rel (p, P)| /m*, and the total 
scattering cross section a{p^^). Denoting the improper 
momentum eigenvectors of tracer and gas by |P) and |p) 
we have 

r = JdPdpng,sv{p,P)a{Te\{p,P))\P){P\®\p){p\. 

Indeed, for separable particle-gas states the expectation 
value of r yields the average total collision rate experi- 
enced by the Brownian particle. 

Let us now see how a single collision changes the 
motional state of the tracer particle according to scat- 
tering theory. If tracer and gas are uncorrelated be- 
fore the collision the outgoing tracer state is p' — 
trgas (S [p (8) pgas] S''^) , where S = I + zT is the two-particle 
scattering operator and a partial trace over the gas has 
to be performed. Employing unitarity, S^S = I, one can 



express the change of the state Ap = p' ~ p as 



Ap = -trgas ([T + Tt,p( 

+ trgas (T [p (g) 



(4) 



"2 trgas (T'^T [p (g) Pgas] + [p<^ Pgas] T^^T) 

The first term generates a constant coherent modifica- 
tion of the unitary evolution and can be absorbed in the 
Hamiltonian H. This energy shift due to "forward scat- 
tering" is usually accounted for by a modified index of 
refraction (l^ and can be disregarded since the gas den- 
sity is uniform. 

The momentum representation of the remaining inco- 
herent part of Q can be expressed entirely in terms of 
the kernel 

(P| trgas (T[|Po)(P[,|® Pgas] Tt)|P'). (5) 

If the momentum-diagonal representation of the gas state 



Pe 



dPoAi(Po) IPoXPoh (6) 



is inserted the kernel (0) assumes the simple form 
(5(P-Po-P' + Po) J(P,P';P-Po), where the 
function 

J{P,P';Q) = /dp„M(Po) (7) 

X (rel (po - Q, P) |To| rel {p^,P-Q)) 
X (rel (po, P' - Q) ITJI rel (po - Q, P')) 

is given in terms of the single particle operator Tq for 
the relative coordinates. Its momentum matrix elements 
are related to the scattering amplitude by (py|To|Pj) — 

{■kK)^^ 5 {p^j ~ pf^ f {pj,p^ 0]. However, inserting 

this into Q one arrives at an ill-defined expression, which 
involves the arbitrary normalization volume 17 and, for 
for P — P' , the square of a i5-function. As is well un- 
derstood, the reason for this is our use of a momentum 
representation of pgas in Plane waves are not in the 
domain of the operator So = lo + *To, which maps in- 
coming asymptotes to outgoing ones (and by extension 
leaves the outgoing ones invariant). 

A possible resort would be therefore to choose a repre- 
sentation of Pgas which permits a decomposition into in- 
and out-states. Indeed, a decomposition into Gaussian 
wave packets admits a conventional but tedious calcula- 
tion of the state change in the limiting case of an infinitely 
massive Brownian particle, m/M — 0, as demonstrated 
in i8j. In the same article it is shown that the identical 
result can be obtained directly by keeping the diagonal 
representation. In this case the extension of Tq beyond 
the physically acceptable domain of incoming states must 
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be complemented by a consistent, physically motivated 
replacement rule: 



(27r;i) I ,2 
^ \{Pf\^o\Pi)\ 



Pt 



\f{Pf^P.)\ 

(^{Pi) \P;\ 



<8) 



In the present case of a finite mass ratio this yields a 
well defined kernel Q for P = P' , and thus implies 
J dQJ {P + Q, P + Q;Q) = 1, which ensures the con- 
servation of the norm of p' . For P ^ P' an extension of 
the rule ^ to different pairs of incoming and outgoing 
relative momenta is required. It can be constructed at 
no additional cost if the momentum change of the two 
pairs is the same, as is the case in 



The temporal change dtp is obtained from Ap in Q by 
replacing T with Tf^/^, 

dtp = trga. (Tri/2 [p® pg,,] fi/^jt) 

-^trga.(ri/2TtTri/2 [p^pgj) 

-^trgas([p®Pgas]rl/2TtTrl/2). (IQ) 

This is suggested by a quantum trajectory unravelling 
of dtp 19] where each trajectory is first weighted by the 
probability of a collision event to take place in an in- 
finitesimal time interval before being scattered. To ob- 
tain the momentum representation of H1U|) we have to 
evaluate 



X 



{2TThY 

n 



{Pf + q\To\Pi + q){p^ - q\Tl\Pf - q)g (q) (P| tr.^s {jr'^' [\Po){Po\ ^ Pgas] rV^jt) |P') (11) 



with arbitrary positive function g. Inserting the formal 
square root of the replacement rule JHJ yields the square 
root of a product of two energy conserving ^-functions 

2 _ 2 

with arguments ± (py — p^) • q. They imply that 

the parallel component q|| = 9||(p _p.) of the momentum 
separation must be zero, which restricts q-integrations to 
the plane perpendicular to the momentum change Pf—Pi- 
This restriction is equally effected by replacing the vec- 
tors q with their projection q^^ = q q|[ onto that plane. 
This way the form of the arguments in each individual 
scattering amplitude already ensures the conservation of 
energy and one is left with a single proper Dirac function 

5 ^ ^^2^' ) ■ Hence, as a natural generalization of 
have 



Using again the momentum diagonal gas 
state ^ reduces the kernel Hll|l to the form 
S{P-Po-P' + P'o)Mi^{P,P';P-Po). It fol- 
lows that the momentum representation of H1U|) is 

dtp (P, P') = JdQ Mi„ (P, P':Q)p{P- Q, P' - Q) 
-\ [Mt, (P) -I- Mt, (P')] p (P, P')(12) 



with 



we 



Mn {P,P';Q) 



{2'Kh) ngas 



X 




f {pf + q±,Pi + q±) 
q^.Pi - q\) 



\J<^ {pi 



\Pi - q±\ 



fgiq±) 



(9) 



It turns into the known replacement rule as q ^ 0. The 
only freedom in this construction is a possible phase fac- 
tor from taking the roots, but symmetry considerations 
lead to the above choice of no additional phase. 

Being able to evaluate traces with a momentum- 
diagonal representation of the gas, we can now combine 
the operator T for the effect of a single collision with 
the rate operator f to obtain the temporal evolution (in 
interaction picture putting aside the contribution of H). 



xV\p^ + q\(^ {Pi + q)^/\Pi - 9k {Pi - q) 

X {pf + q|To|p, + q){p, - q\ll\Pf - q) 

and M„=it (P) := / dQ A/i„ {P + Q,P + Q;Q). Here I 
introduced Pj :— rel (^Pq, — and Pj ■— Pi ~ Q 

as functions of Pq, and defined q :— rel ^0, ^ . 

In order to evaluate Mjn the integration is now trans- 
formed to dpj. Incidentally, this suggests a natural fac- 
torization of the p distribution into a product of square 
roots, ^J'J^(p^^)J^{p^ , since Pq can be equally expressed 
as a function of P or of P'. Applying the replacement 
rule © projects q to = qj^(^p^_p -j not only in the 
scattering amplitudes, but also in the argument of p. 
One obtains the well-defined expression 
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M ^ m* / V 7\/ in 



As the last step, the transformation p^ K 



— — — Y factorizes the integrand into P and 



M 2 m, 2 

P' contributions, 
Mi„ (P, P'; Q) = y dii:^ {K-Q)F {K, P-Q;Q) 
xF*{K,P'-Q;Q) (13) 

with given in By returning to the Schrodinger 
picture and noting / dK6 {K ■ Q) = Q^^ Jqj_ dK one 
finds that p2|l with (|13l) is the momentum representation 
of (^), which closes its derivation. Note that the form of £ 
is a consequence only of the premises and the replacement 

The physics described by ® is easy to discuss in the 
momentum basis (|12ll . We will see that M^^j (P) is the 
rate of a classical particle with momentum P to be scat- 
tered by the gas into a different direction or velocity. 
Hence, the second line in H12|l effects a reduction of the 
coherences p (P, P' ^ P) determined by the arithmetic 
mean of the corresponding momenta on the diagonal. 
The first line in (|12|l . on the other hand, may reduce 
or produce coherences. While the complex quantity 
cannot be decomposed into classical rates, it is reassuring 
that it involves an integration over all possible, in gen- 
eral nonparallel, pairs of two-particle scattering trajecto- 
ries which end at the Brownian momenta {P,P'), each 
part conserving the total energy and the momentum with 
an exchange of Q, and weighted by the thus restricted 
distribution of available gas momenta. [Diosi's equation 
|ll|. in comparison, involves the differential cross section 
da/d^l = l/l^ and therefore at most pairs of identical 
trajectories.] 

After a time long compared to the time scale of de- 
coherence the motional state is expected to be practi- 
cally indistinguishable from a classical state. As such it 
should be characterized by the momentum distribution 
w (P) = p (P, P) alone, and one expects that the motion 
of the diagonal elements predicted by is equal to the 
classical linear Boltzmann equation. Indeed, one obtains 
from (112) 

dtw{P) ^ j dQMi{P,Q)w{P-Q) 



where Mf^^ (P) = / dQM^^ (P + Q; Q) and 



da , 



(15) 



is the rate density of a classical Brownian particle to be 
scattered into momentum P upon a momentum exchange 
of Q. Here p^, = rel [K, P Q), p^j = p^, - Q, and 

da/dn = l/l^. Note that the classic form of (dj) is 
obtained by transforming the Q-integrals over the Dirac 
5 into angular integrations over the relative momentum 
direction. It follows that the stationary solution of 
is given, for thermal gas states /x, by the correspond- 
ing (momentum diagonal) thermal Brownian state of the 
classical equation, and that the H-theorem applies. 

Another border case of is the limit of an infinitely 
massive tracer. A short calculation confirms that letting 
m/M approach zero reduces Lp to the (corrected ver- 
sion of the) master equation by Gallis and Fleming 
, which attributes the loss of coherence to the amount 
of position information gained by the colliding gas. Fi- 
nally, assumes the form of Vacchini's equation 
if one replaces the true scattering amplitude / in Q by 
its Born approximation /b (which depends only on the 
momentum transfer). 

In conclusion, I presented the quantum version of the 
linear Boltzmann equation. It unifies, in the form of a 
completely positive master equation, the decohering and 
dissipative dynamics in the motion of a Brownian parti- 
cle, and it comprises various known dynamic behaviors 
as limiting forms. 
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-Mt, {P)w{P)., 



(14) 
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